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Abstract 

We study saturation effects in the production of dijets in p-p and p-Pb collisions using 
the framework of high energy factorization. We focus on central-forward jet configurations, 
which allow for probing gluon density at low longitudinal momentum fraction. We find 
significant suppression of the central-forward jet decorrelations in p-Pb compared to p-p, 
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' which we attribute to saturation of gluon density in nuclei. 

m ' 

O 

in ; 1 Introduction 

in; 

Physics in the forward region at the Large Hadron Collider (LHC) is a very interesting and 
exciting field since it involves interplay of the kinematical scales, like for example transverse 
momenta of jets, with scales generated by the QCD dynamics, like the saturation scale Q s [T]. 
The latter scale characterizes formation of dense system of partons and there is a growing 
evidence that the phenomenon of saturation of gluons indeed occurs [2H5]. To further advance 
studies of saturation and other possible effects occurring at high partonic density, the LHC is 
going to collide p-Pb this year. This will allow, in particular, for the study of the onset of 
saturation as a function of variables related to the transverse momentum pt of jets. In addition, 
the understanding of the interplay of scales in the jet production in p-p and p-Pb will permit 
to constrain the unintegrated parton densities in the large phase space available for partons, 
i.e. 10~ 6 < x < 0.1, 5GeV < kt < 150 GeV, where x is the longitudinal momentum fraction of 
the hadron carried by a parton, while kt is the component of its momentum transverse to the 
collision axis. In particular, the study of exclusive final states, like jets, allows for determination 
of the unintegrated gluon density in the range of large momenta. 

There are formalisms that allow one to study dense systems [6j or systems with hard momen- 
tum scale involved [7]. However, the formalism which accounts for both the high energy scale 
and the hard momentum scale pt is provided only by the high energy factorization [8]. In this 
framework, one of the elements that enter the factorization formula is the unintegrated gluon 
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density. Depending on approximation, it satisfies the BFKL, BK, or CCFM evolution equa- 
tions [9l416|. The BFKL and BK equations are known already at NLO and sum up emissions 
of gluons with strong ordering in the longitudinal momentum fractions of subsequently emitted 
gluons. The important issue for phenomenological applications, in particular for exclusive final 
states at large pt, is to perform resummations of most relevant higher order corrections to the 
evolution kernel of the BFKL equation. This is because only then, the solution of the equation 
for the unintegrated gluon density is physically relevant and well defined |17H19j . 

Another issue is that since the BFKL or CCFM equation are linear they predict strong rise of 
gluon density at small values of gluon's kt which leads to conflict with unitarity bounds. Effects 
of higher orders, although suppress the growth of gluon density, preserve its power like behavior 
as a function of kt at low kt values. To restore unitarity one supplements the BFKL equation 
with a nonlinear term which accounts for fusions of gluons. These unitarity corrections, which 
are taken into account in the BK or JIMWLK equations, give rise to an emergent semi-hard 
scale, called the saturation scale Q s (x), at which the gluon density has a maximum and which 
therefore defines the most probable momentum of gluon |20p 21j. 

An interesting process in which both saturation and production at high pt can be studied, 
is the process of dijet production where the jets are separated by large rapidity |22H25| . More 
specifically, we shall focus on the case in which one jet is in the central while the other in 
the forward rapidity region. Such a final state probes parton density of one of the protons at 
low longitudinal momentum fraction x while the other at large longitudinal momentum fraction. 
The latter proton can be described by the well known collinear parton distribution functions and 
therefore such process is perfectly suited to study properties of the unintegrated gluon density 
at low x and especially its saturation. 

In this paper, we discuss production of dijets in the p-p and p-Pb collisions. The former 
serves as a benchmark, i.e. first we fit the unintegrated gluon density to the F 2 data [26j and then 
apply it to calculate observables in p-p characterizing the dijet system like angular correlations 
of produced jets, pt spectra of forward and central jets and their rapidity distributions. In the 
next step we compute predictions for rapidity distributions and angular decorrelations of central- 
forward jets produced in p-Pb collision. For both observables we see significant suppression of 
the cross section due to saturation effects in the nucleus. 

The paper is organized as follows. In section [2] we introduce the high energy factorization 
framework and define the observables we want to use as a tool to study saturation effects. In 
section [3] we introduce the unified BK/DGLAP equation |27| for the unintegrated gluon density 
and present results of fits of the unintegrated gluon to the combined HERA data. In section [J] 
we apply the high energy factorization framework together with our fitted unintegrated gluon 
to calculate observables for central-forward dijet system. In the section [5] we calculate rapidity 
distribution and angular decorrelations of central-forward jets produced in p-Pb collision. We 
conclude our studies in section [6j 

2 Dijet production in high energy factorization approach 

The main goal of this paper is to provide predictions for azimuthal decorrelations of jets produced 
in p-p and p-Pb collisions. Consider the process of the production of a dijet system in the collision 
of two hadrons 



A + B -> Ji + J 2 + X . 




The leading order contribution comes from the 2 — > 2 partonic process 
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Figure 1: Jet production in the forward region in hadron-hadron collisions. 



In this study we focus on the asymmetric configuration with one jet produced in the forward and 
the other in the central rapidity region as illustrated in Fig. [TJ The fractions of the longitudinal 
momenta of the initial state partons are related to the transverse momenta and rapidities of the 
final state partons by 

xi = -±= {p tl e^ + p t2 e y >) , x 2 = -±= (ptie- y ' + Pt2 e~ y2 ) , (2.3) 

where S is the squared energy in the center of mass system of the incoming hadrons. Hence, our 
central-forward configuration corresponds to one of the XjS being small and the other one large. 

The first study of such configurations in the p-p collisions was performed in [28H30] using 
the CASCADE Monte Carlo generator |31| . Assuming, without loss of generality, that x\ ~ 1 
and x 2 <C 1 we have 

K = x i p i, ( 2 - 4 ) 
k% = x 2 P2 + k?, (2.5) 

where we used the Sudakov decomposition of the initial partons' 4-momenta. Here P^ 2 are 
the 4-momenta of incoming hadrons, which, in the center of mass frame, take the form P^ 2 = 
(1, 0, 0, ±1) and P\-P 2 = ^S. The momentum of the off-shell parton satisfies k 2 = k 2 = —k 2 , 
where k = kt = |k|. This leads to the following form of the cross section 

, , ^ , 2 = Y] 1fi 37 1 ^ 2 M a g^cdXlfa/A{xi,H 2 )(/) g /B(x 2 ,k 2 ,H 2 ) ^ , (2.6) 

dyidy 2 d 2 p u d z p 2t ^ d lQ-K 4 {xix 2 Sy a 1 y/ 1 + 5 cd 

and 

k 2 = Pti + Pt2 + 2 PtiPt2 cos Acj) , (2.7) 

where A<p = (fti — 4>2 is the azimuthal distance between the outgoing partons and M. ag _^ C( i is 
the matrix element for the 2—7-2 process with one off-shell initial state gluon and three on-shell 
partons, a,c,d, which can be either quarks or gluons |32j . The following partonic sub-processes 
contribute to the production of our dijet system 

qg^qg, gg^-qq, gg^gg- (2-8) 

On the side of the off-shell gluon in Eq. (|2.6p . we have the unintegrated gluon density </> 9 /_b(x2, k 2 fx 2 ), 
which depends on the longitudinal momentum fraction x 2 , on the transverse momentum of the 
off-shell gluon, and in general as well as on hard scale fi. The hard scale dependence introduces 
DGLAP-like ordering effects in the high energy factorization framework and makes it applica- 
ble in studies of exclusive final states. In our calculations, however, we follow the KMS [UJ 
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Figure 2: Diagrammatic representation of the F 2 structure function in high energy factorization. 

scheme to introduce corrections to the gluon density which make it applicable to the studies of 
jet physics. Because of this we shall skip the argument [i in the expressions for the unintegrated 
gluon density below. 

On the side of the on-shell parton, which is probed at high values of the longitudinal mo- 
mentum fraction x\, it is legitimate to use the collinear parton density f a /A{ x i^ 2 )- 

The above result depends only on the difference of the azimuthal angles A<j), so one can 
change variables and integrate out one of angles (pi. This leads to 

, , , d(T , Ta~T = B 2< tim G ^2 M ^cdXlfa/A^l, ^) <Pg/B{x 2 ,k 2 ) 1 , (2.9) 

dyidy 2 dpi t dp 2 tdA(t) ^ 8ir z (xix 2 Sy a 1 y > 1 + 5 cd 

with k 2 = p\ x + p 2 2 + 2ptip t2 cos Acf>. 

3 Unintegrated gluon density from the unified BK/DGLAP frame- 
work fitted to combined HERA data 

The formulation of the NLO BFKL equation [33— 35J has been known already for some time. 
Also the NLO BK equation has been derived |36| but, because of its complicated structure, 
only solutions of some approximate forms of the BK equations are known (see |37H40| ). The 
basic formulation of the NLO BFKL equation is unstable (due to non-positive definite kernel) 
and in order to stabilize it one needs to resume a subset of higher order corrections |17 | ll8 | l4"T]. 
In our study, we will use the approach to this problem formulated in |41| in which large part 
of the higher order corrections is provided by the consistency constraint on emissions of real 
gluons. The other important corrections are coming from running of the coupling constant 
and the nonsingular pieces of the DGLAP splitting functions. Other approaches were discusses 
in [HH]. 
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Figure 3: The proton structure function F2(x,Q 2 ) from the fit of our framework, in its linear 
and nonlinear variant, to the combined data from HERA |26| as a function of x for the Q 2 range 
from 1.5 to 400 GeV 2 (with the vertical offsets of 0.2). 



The corresponding equation for the unintegrated gluon density reads |27y45| 

<P p (x,k 2 ) = cPP(x,k 2 ) 

a s (k 2 )N c [Uz r dl 2 { Z 2 P (f,Z 2 )fl(f -I 2 ) - k 2 4> p (l,k 2 ) | k 2 ^k 2 ) 

|4Z 4 + fc 4 |t 
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2aj(k 2 ) 
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P " {Z) " 2 ~ir) (~ z ,l 2 ) + zP 9q (z)V 

°° dl 2 , , , 2 \ 2 2 , f°° dl 2 , fl 2 \ 



k' 2 



I 2 



4> p (x,l ) + 4> p (x,k ) 



k 2 



p ln ( ) M x > 1 ) 



(3.1) 



where z = x/x' (see Fig. [2] for explanation of the variables). For convenience, we omit the g 
subscript in the unintegrated gluon density symbol and keep only the subscript denoting the 
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Figure 4: Left: Unintegrated gluon density in the proton from the solution of Eq. (I3.ip with free 
parameters fitted to F2 data. Middle: Comparison of the gluon density from the left plot with 
the unintegrated gluon from the fit using only the linear part of Eq. (|3,ip . Right: Unintegrated 
gluon density in the Pb nucleus obtained as for the proton case but with the nonlinear term in 
Eq. enhanced by A 1 / 3 with A = 207. 

hadron. The theta function in Eq. (|3.ip introduces the kinematical constraint and the nonlinear 
term, which supplies unitarity corrections, is given by the triple pomeron vertex |46| . The 
two terms in the third line in Eq. (|3.ip correspond to the DGLAP effects generated by that 
part of the splitting function P gg (z) which is nonsingular in the limit z — > and by the quarks 
respectively, with k 2 ) corresponding to the singlet quark distributions. The At LO in In 1/x 
this equation reduces to the BK equation, after performing Fourier transform to the coordinate 
space |13 [[2Tll4'6ll47] . For numerical method to solve Eq. (|3.ip we refer the reader to |27j . 

The strength of the nonlinear term in Eq. (|3.ip is controlled by the parameter R which has 
an interpretation of the proton radius and it comes from integration of the gluon density over 
the impact parameter b. In our framework, we assume the uniform distributions of gluons in 
the nucleon therefore our gluon density is proportional to Q(R — b) |27j . 

The input gluon distribution <p p °\x, k 2 ) is is given by 

where xg(x, fcg) is the integrated gluon distribution at the initial scale, which we set to k 2 = 
lGeV 2 , and we take the following parametrization 

xg(x) = N{1- xf{l- Dx) , (3.3) 

which is similar to what was used in [48J. For the strong coupling, we take the one- loop result 
with Aqcd set to 350 MeV. 

The evolution equation (|3.ip is used to determine the unintegrated gluon above the initial 
momentum scale that is for k 2 > 1 GeV 2 . In the region k 2 < 1 GeV 2 , the gluon density (j) p (x, k 2 ) 
is constrained by the condition that it should match the evolved unintegrated gluon density at 
k 2 = k 2 . We choose to parametrize the unintegrated gluon in this region by 

<p p {x, k 2 ) = k 2 <p p (x, k 2 = 1 GeV 2 ) for k 2 < 1 GeV 2 , (3.4) 

which is motivated by the shape of the gluon density obtained from solution of the LO BK 
equation in the saturated regime |49j. 

The unintegrated gluon density from Eqs. (|3. 1|) (|3.3|) convoluted with impact factors taken 
from [H] allows one to compute the structure function F%(x, Q 2 ), which in turn can be used to 
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fit the free parameters of the model. We performed such a fit to the combined HERA data |26| 
in the kinematical range of x < 0.01 and the full range of Q 2 . As shown in Fig. [3] (red solid 
line), we obtain a very good description of data, which corresponds to x 2 /ndof =1.73 and the 
following values of the parameters: N = 0.994, j3 = 18.6, D = -82.1 and R = 2.40 GeV -1 . 

In Fig. 2] (left) we show the unintegrated gluon density, corresponding to the above fit, as a 
function of the transverse momentum of the gluon for a range of x values. The gluon from the 
evolution is sewed at k 2 = 1 GeV 2 with the parametrization (|3.4p . The sharp peek corresponds 
to the point in k where the matching was done. We see however that as one goes to lower x 
values, perturbatively generated maximum starts to emerge. That is a signal of the presence of 
saturation scale defining the most probable momentum of the gluon. The results for estimated 
value of the saturation scale are in an agreement with |40 | l4*9" l l50] 

Our main focus in this study is on the framework with saturation of gluon density described 
above which is based on the nonlinear evolution equation (|3.ip . It is however interesting to 
compare our results to the case in which the gluon is determined from the framework without 
saturation. This is naturally provided by the linear version of Eq. (|3.ip that corresponds to 
dropping the last term on the right hand side of that equation, which now becomes independent 
of R. We performed an analogous fit to the one described above but taking the linearized version 
of Eq. (|3,ip and restricting Q 2 to the values above 4.5 GeV 2 to stay outside of the region where 
the saturation effects may be important. The fit parameters at the minimal value of x 2 = 1-51 
are: N = 0.004, j3 = 26.7 and D = —51102. The corresponding results are shown in Fig. [3] (blue 
dashed line) for the whole range of Q 2 including the bins below 4.5 GeV 2 , which were not used 
in the fit. We see that the linear gluon gives too strong rise of F2 with x especially at low values 
of Q 2 . This remains true even if we fit the linear version of Eq. (|3.1|) to the full range of Q 2 . 
The best value of \ 2 we were able to achieve in this case was 3.86. We therefore conclude that 
some mechanism damping the gluon density at low x and low Q 2 is necessary to describe the F2 
HERA data in the full range of Q 2 . This mechanism corresponds to saturation of gluon density. 

We summarize this part of our study by comparing the two versions of the unintegrated gluon 
density, linear and non-linear, in Fig. [?] (middle) . The linear gluon for k 2 > 1 GeV 2 corresponds 
to the fit described above. For the values of k 2 below 1 GeV 2 , similarly to what we did for the 
non-linear case, we parametrize our gluon, this time by <p(x,k 2 ) = <p(x, fcg = lGeV 2 ). We see 
in Fig. 3] (middle) that at large values of x and kt, both distributions are similar. As one goes 
to lower kt, however, the linear gluon rises much faster than the non- linear one. This effect 
becomes significantly stronger for smaller values of x. 

4 Central-forward dijet production in p-p collisions at the LHC 

We are now in the position to compute the cross section for central-forward jet production in the 
proton-proton collision. We use directly the formula (|2,9p with the matrix element calculated 
in |32j and the unintegrated gluon distributions determined in the previous section. For the 
collinear parton distributions we take CTEQ6mE |51j. 

The CMS collaboration has reported the measurement |52j of the transverse momentum 
distributions of pairs of jets, one of which is restricted to be in the forward and the other in the 
central part of the detector. This measurement has already been confronted with various Monte 
Carlo predictions including HE J |53t l54"] and CASCADE |31j. We compute the corresponding 
distributions in our framework taking the selection cuts which match those of CMS. We set 
the normalization and factorization scale equal to the fixed value of 60 GeV and vary it by the 
factors 1/2 and 2 to asses the uncertainty. 

The results for jet pt spectra are shown in Fig. [5] for the central (left) and the forward (right) 
jet respectively. We find that our predictions based on nonlinear (red) and linear (green) evolu- 
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Figure 5: The transverse momentum spectra of the central (left) and forward (right) jets 
from our framework compared to the CMS data. The uncertainty bands come from varying the 
factorization and renormalization scale by factor two around the central value of /x = 60 GeV. 



tion equations reproduce the pattern of the CMS measurement. As expected, this observable in 
experimentally accessible region is weakly sensitive to saturation. We see some suppression of 
the cross section based on the nonlinear evolution equation which is a manifestation of the lower 
gluon density in the small kt region. The description with both densities is good despite the fact 
that our modeling of jets is very simple since, in the framework we use here, each of the two 
jets is just a single parton. Moreover, the small excess of low-p^ is expected precisely because 
of that. Adding a parton shower on the top of our partonic result would cause some energy to 
go out of the jet and that would affect the low-pt jets since they are broader then the high-p^ 
jets. For related discussion within the HE J approach see [54^55] . Similarly, we may expect 
that due to large rapidity gap between the produced jet, some extra BFKL-type radiation could 
be emitted. An explicit demonstration of the role of these effects in our framework opens an 
interesting possibility for future work. However, they are not crucial for the following discussion 
since they will affect the results for p-p and p-Pb in the same way and our focus in this work is 
on the relative differences between this two cases. 

Also the slightly smaller cross section in the high-pt region for the central jet is well un- 
derstood. As follows from Eq. (j2.3|) . the value of xi probed by the central jet pt ~ 140 GeV 
corresponds to X2 ~ 0.02 which is beyond the limit we used in our fits of unintegrated gluon 
and therefore our predictive power in this region is limited. We have checked that if the gluon 
is fitted with the upper limit on x extended to 0.02, the result for pt distribution of the central 
jet becomes consistent with the data also in the high-pj region. Trying to extend our framework 
to higher values of x would be in itself an interesting project. It goes however beyond the scope 
of this paper, therefore throughout this study, we restrict ourselves to the region of x < 0.01. It 
should be emphasized, that even without going to large values of x the distributions from Fig. [5] 
are unique since they give direct access to the unintegrated gluon at medium and large values 
of kt, provided that one works in the framework which permits to compute gluon in that region. 
The framework we adopted for this study satisfies that criterion. That opens the possibility 
to study the importance of terms of higher order from the point of view of BFKL, i.e. energy 
conservation and subleading pieces of the splitting function. 
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5 Signatures of saturation in central-forward dijet production in 
p-Pb collisions 

In the preceding sections, we have shown that the high energy factorization formalism with the 
unintegrated gluon from the QCD evolution equation with saturation can successfully account 
for the features measured both in the e-p and p-p collisions. The slightly less precise description 
of the data in the latter case is expected and fully understood, and it can be traced back to the 
purely partonic nature of our result or to the upper limit on the x used in our fits of unintegrated 
gluon. 

One of the main features of the formalism we use is that the unintegrated gluon is determined 
from the nonlinear evolution equation and therefore it exhibits saturation effects around certain 
momentum scale Q s (x). In this section, we address the question whether those effects could be 
studied experimentally in collisions at the LHC. 

The main challenge of such a study, as pointed out by numerous analysis of HERA data |4CH 
I56l - I59| . is that the saturation scale in the proton seems to be of the order of a few GeV, hence it 
lies just at the border between perturbative and non-perturbative regime of QCD. This makes 
it difficult to access both theoretically and experimentally. One way to improve the situation 
is to go to the p-A collisions since, as widely discussed in the literature, the saturation scale in 
the nucleus is expected to be significantly higher than in the proton [31ElEDrE2]- To estimate 
the possible effects of saturation in the heavy nucleus we use a simple formula for the nucleus 
radius following from counting the number of nucleons for Woods-Saxon nuclear density profile. 
The radius of the nucleus reads 

R A = RA 1/3 , (5.1) 

where R is the proton radius, which is one of the fitted parameters of our framework as described 
in section [3] and A is the mass number (A = 207 for Pb, A = 196 for Au). The above definition 
has the property that in the limit A — > 1 the result for the proton is recovered. 

Analogous equation to Eq. (|3.ip for the Heavy Ion (HI) normalized to the proton reads 

4>m/A(x,k 2 ) = (j)^l /A (x,k 2 ) 

a s (k 2 )N c f 1 dz f°° all 2 / l 2 0m/A^,l 2 )O(^-l 2 ) ~ k 2 <p m/A {%k 2 ) k 2 <f> m/A ^,k 2 ) 



7T J x Z 

2\ /■! 



a s (k 



kl I 2 I \l 2 -k 2 \ |4Z 4 + fc 4 |5 

k 2 



27Tk 2 

2A l l*a 2 s {k 2 ) 
R? 



dz (p gg ( z )-^j dl 2 fa yA (~J 2 )+zP gq (z)X HI/A (~,k 
I J 2 "prill/ A fa 1 ) 1 +(f>m/A{x,k ) J ^ -p In I ^ I 4>m/A\ x -> ' ) 



(5.2) 



where we used Eq. (|5.ip to express the radius of the heavy ion in terms of the proton radius 
thus (f>ui(x, k 2 ) = A<f>Y{i/ A {x, k 2 ) is the distribution of gluons per nucleon in the nucleus. We see 
that the strength of the nonlinear term in Eq. (|5.2[) is enhanced by A 1 ^ 3 (6j,[63]. We are aware 
that this modification is not sufficient to fully model the nuclear target and is a somewhat crude 
approximations but it will suffice as a first approximation to estimate the saturation effects in 
the nucleus. For a more sophisticated approach see |64| and references therein. 

In Fig. [|] (right) we show the gluon density in the Pb nucleus which results from the ap- 
plication of the above prescription to the unintegrated gluon distribution in the proton from 
Fig. [3] (left). We notice that due to stronger saturation effects in Pb, the gluon density is lower 
than that in the proton for the same value of x. We also see that in Pb the maxima are shifted 
towards larger values of the gluon's transverse momentum kt, which corresponds to the larger 
saturation scale in the nucleus. 
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Figure 6: Ratio of d-Au/p-p coincidence probabilities C(A<p) for the forward dihadron pro- 
duction at RHIC as a function of the azimuthal distance between the particles. The d-Au data 
were shifted by a constant. The errors of the ratio were determined from relative errors of each 
C(A(p) before the shift. The band corresponds to our prediction with the uncertainty related 
the unknown yield of uncorrelated dihadron production. 



An observable which is very well suited to study saturation is the azimuthal correlation of 
the central and forward jet. It is an inclusive observable that measures radiation between jets 
and is therefore sensitive to potential saturation effects which are supposed to decrease the rate 
of emissions when the parton density is probed at very low x where the high density of partons 
leads to their recombinations. The most interesting region of Acj) is that close to ir since the 
produced jets are almost back to back and the gluon density is probed at low kt- 

Before we turn to the discussion of our main results for the central- forward dijets production 
at the LHC, it is interesting to check how the magnitude of the suppression of decorrelations 
predicted by our framework compares to the existing data for the forward dihadron production 
from RHIC. STAR and PHENIX measured the coincidence probability defined as the ratio of the 
yield of the 7T°-pair to the inclusive 7r° yield, C(Acf)) = A pair (A(/>)/7V incl . Precise determination 
of this observable requires convolution of our diparton pair from Eq. (|2,9p with the fragmentation 
functions as well as computation of inclusive ir° production and also the uncorrelated tt° pair 
production (see [3,61j for corresponding results in different frameworks). All this goes beyond 
the scope of the present paper and will be published elsewhere. Nevertheless, we can make a 
meaningful estimate for a related quantity that can be compared to the STAR data [651166} . 
To minimise the contribution from the uncorrelated pair production we shift the d+Au data 
by -0.00145, similarly to what has been done in [3], so that the d+Au and p-p data coincide 
below A(f) = tt/2. Assuming that the remaining yield comes predominantly from the correlated 
dihadron production, we compare the ratio C dAu (Aiji>)/C pp (A(/>) for the experimental data with 
that from our computation for diparton production at ^/s = 200 GeV with pti > 2.5 GeV, 
Pti > Pt2 > 1-5 GeV and 2.4 < y\i < 4.0. By using the above ratio we do not need to worry 
about the normalization to the inclusive spectra and it is also reasonable to expect that most of 
the effects from parton fragmentation will cancel. Our absolute prediction (with = pti) for the 
ratio of the coincidence probabilities is shown in Fig. [6] together with the data from STAR [6 5"|I66| . 
Since our framework does not depend on the impact parameter, we chose to compare to the 
data averaged over centralities. The band in Fig. [6] corresponds to the uncertainty related to the 
assumption of negligible contribution form the uncorrelated production after the shift described 
above. To asses this uncertainty we considered an alternative scenario in which half of the yield 
seen in the data at A<p = tt/2 is attributed to the independent production and the other half 
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Figure 7: Differential cross sections for central- forward dijet production as functions of az- 
imuthal distance between the jets Acfi (top) or rapidities of the jets (bottom) for the case of p-p 
and p-Pb collisions and three different cuts on jets' pt- 



to the correlated production that we can predict with our framework. As we see in Fig. [6j the 
suppression pattern of the away-side peak of the dihadron spectra from d-Au collisions at RHIC 
is correctly reproduced by our calculation which shows that our theoretical framework captures 
the essential physics of this class of processes. 

We move now to the central-forward dijet production in the p-Pb collisions at the LHC. In the 
top row of Fig. [7] we show the differential cross section for the central-forward dijet production as 
a function of the azimuthal distance between the jets for the p-p and p-Pb collisions. To obtain 
those results we employed the linear and nonlinear versions of the evolution equation (|3,ip for 
the proton and the evolution equation (|5.2p for Pb. We used selection similar to that form the 
previous section except for the pt cut which we now vary from 15, through 25, to 35 GeV and the 
rapidity which is restricted to positive values. The latter corresponds to the fact that, contrary 
to the p-p case, the p-Pb collision is asymmetric and, as follows from Eq. (|2,3p . one probes the 
gluon density in Pb at low x only by measuring the forward jets going into the region of positive 
rapidity. 

First observation from Fig. [7] is that the non-linear evolution leads to a significant suppres- 
sion of the Acj) and rapidity distributions already for the proton case. This alone is a clear 
manifestation of saturation. Then we see that the A(f> cross section near the peak region is 
suppressed further by the factor of about two for the case of the p-Pb collision and the effect 
extends to lower values of Acft as we lower the pt threshold (going from right to left plot). This 
is precisely the consequence of gluon saturation which is stronger in the Pb nucleus then in the 
proton and therefore the unintegrated gluon distribution in the region of small and medium kt 
is suppressed in Pb compared to the proton case as shown in Fig. 0] (right). It is this region of 
gluon's kt that is probed by the dijet configurations with A(p ~ tt and that is what leads to the 
lower cross section in the area of the peak. 
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Figure 8: Differential cross sections for central- forward dijet production at yfs = 5 TeV and 
8.8 TeV as functions of azimuthal distance between the jets Atp for the case of p-p and p-Pb 
collisions and three different cuts on jets' pt- 

We notice that the non-linear results have a dip structure near Acf) ~ ir. This is a consequence 
of the feature of a high energy factorisable gluon density which goes down to zero like k 2 . On 
the other hand, as discussed in section [3J in the linear case we model the behaviour of the 
unintegrated gluon density by assuming that the gluon density behaves like a constant and 
therefore the linear result for the Acp distribution keeps growing as Acf) — > ir. These features make 
the Acf> distribution a particularly interesting observable for testing shapes of gluon densities and 
more generally the validity of the high energy factorization as pointed out recently in [67] , On 
the top of that, the region near Acp ~ ir is also sensitive to Sudakov (virtual corrections) and 
parton shower effects (taking energy from the jets we measure) which have the tendency to 
reduce the cross section in the region Acp near the peak. Some refinement along those lines 
could be envisaged in the future. These effects act however in a similar way for the proton and 
for the heavy ion since they affect the hard scattering. Because we are interested in searching 
for saturation effects in the initial state parton density, it is legitimate to neglect them in this 
study and focus on the relative difference between the cases of p-p and p-Pb collisions. The 
main point we would like to emphasize here is that the suppression due to saturation predicted 
in Fig. [7J is both strong and it extends for large enough range in Acf) to allow for experimental 
discrimination between the linear and non-linear scenario, even if the very small region near 
Acp = 7r may profit from further refinements. 

In the bottom row of Fig. [7] we present the rapidity distributions of forward and central jets 
for the case of p-p and p-Pb collisions. As expected, the saturation effects which are stronger 
in the nucleus than in the proton lead to lower yields both for the central and forward jet 
production in the p-Pb collision. Consistently with the decorrelation results, also here, the 
difference between p-p and p-Pb becomes more pronounced as one lowers the value of the jet pt 
cut. 

Finally, in Fig.[8]we show decorrelations plots similar to those from Fig. [JJbut for the energies 
of the actual p-Pb collisions, i.e. the current y/s = 5 TeV and the nominal \fs = 8.8 TeV. As 
expected, the total yields increase with energy but the relative difference between the p-p and 
p-Pb case seems to remain similar. 

6 Conclusion and outlook 

We presented the analysis of e-p, p-p and p-Pb collisions in the framework of high energy 
factorisation with the unintegrated gluon density given by the nonlinear QCD evolution equation. 
We have shown that this formalism can successfully account for features measured in e-p and 
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p-p data. For comparison, we also performed calculations within the linear evolution framework 
and discussed differences between results from the two scenarios. 

We then used the above non-linear framework to provide an estimate of the effects of gluon 
saturation in the nuclei. We presented predictions for the azimuthal decorrelations as well as the 
rapidity distributions for the p-p and p-Pb collisions. Our main finding is that saturation in the 
Pb nucleus has a potential to manifest itself as a factor two suppression of the central-forward 
jet decorrelation in the region of the azimuthal distance between the jets A(J) ~ ir. The effect 
becomes more pronounced with lower cuts on jets' transverse momenta. 

The framework used in our study allows for a number of refinements that would lead to 
a better description of data as well as for more accurate predictions. We could, for example, 
extend our analysis by introducing non-trivial impact parameter dependence of the unintegrated 
gluon density. That would allow us to study saturation effects as a function of centrality of the 
p-Pb collision. 

Another interesting possibility for future work is opened by an ongoing discussion on break- 
down of the high energy factorization and related issue of multiple definitions of the unintegrated 
gluon density. As advocated in [681169] , this generalized description provides a framework which 
is better theoretically motivated. There are two reasons why the study in such a framework 
would be interesting. First of all, because that framework was derived only with simplified 
matrix elements, whereas in our study we have the exact ones. That would allow one to inves- 
tigate the relative importance of the corrections to the high energy factorization formula from 
section [2] compared to the case with exact kinematics. The second reason is that the consensus 
as to which gluon distribution should be used to study dijet production has not been reached 
yet. The definition of the unintegrated gluon that we used in this study follows directly from 
Feynman diagrams and, as argued in |70j, is therefore the valid form of unintegrated gluon den- 
sity. Hence, we performed our study in the framework of the high energy factorization and the 
results from sections U] and [5] can then be used in the future as a benchmark for further studies 
within extended formalisms. 

Finally, we note that the approach used in this study is unique as it allows one to study 
hard final states in the framework with saturation. This offers the possibility to constrain the 
unintegrated gluon density in the large range of momentum available at the LHC. Our current 
limit of x < 0.01 could be extended to higher values of x by using the gluon density from the 
CCFM equation with saturation |71j . That would allow for a better description of the high pt 
tail of jet spectra. 
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